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In this paper we calculated the n-point hard-thermal- loop (HTL) vertex functions in QCD/QED 
for n= 2, 3 and 4 in the physical representation in the real time formalism (RTF). The result showed 
that the n-point HTL vertex functions can be classified into two groups, a) those with odd numbers 
■ of external retarded indices, and b) the others with even numbers of external retarded indices. The 

' n-point HTL vertex functions with one retarded index, which obviously belong to the first group a), 

, are nothing but the HTL vertex functions that appear in the imaginary time formalism (ITF), and 

vise versa. All the HTL vertex functions belonging to the first group a) are of 0(5^r^) , and satisfy 
among them the simple QED-type Ward-Takahashi identities, as in the ITF. Those vertex functions 
^ belonging to the second group b) never appear in the ITF, namely their existence is characteristic of 

. the RTF, and their HTL's have the high temperature behavior of O(g^T^), one-power of T higher 

' than usual. Despite this difference we could verify that those HTL vertex functions belonging to 

I the second group b) also satisfy among themselves the QED-type Ward-Takahashi identities, thus 

guaranteeing the gauge invariance of the HTL's in the real time thermal QCD/QED. 

^ ■ PACS numbers: ll.lO.Wx, 11.15.Tk, 12.38.Cy 

in ■ 
^ , 

INTRODUCTION AND SUMMARY 

The hard-thermal-loop (HTL) resummed effective perturbation theory of Braaten and Pisarski [|l|, ^ has given 



> 

o 



' us in finite temperature (or thermal) field theory ^ a procedure correctly estimating/extracting the dominant 
temperature effect due to the semi-classical thermal fluctuation. HTL's are originally determined in the imaginary time 
O ^ formulation of thermal field theory [Q, |^ by computing the one-loop diagrams. For calculating the Feynman diagrams 
in thermal field theory there are essentially two different formulations, namely the imaginary time formalism (ITF) and 
the real time formalisn(RTF) ||, |^, The ITF is restricted to calculating the static quantities in equilibrium, while 
^ [ the RTF is indispensable for analysing the dynamical quantities and for investigating the out-of-(or non)-equilibrium 
• , system thus for studying the physics of quark-gluon plasma. 

' The relation between quantities calculated in the ITF and in the RTF has been studied so far [^; the n-point 
fimction in the ITF corresponds to a specific sum of the n-point functions in the RTF, the same is true for the n-point 
HTL functions. In the ITF all the n-point HTL functions have been obtained, and are shown to be ultraviolet finite, 
gauge invariant and to satisfy simple Ward-Takahashi identities |^. Also in the RTF the n-point HTL functions 
{n—2, 3, 4) and their spectral functions have been partly determined 0, ^, but the explicit expressions of the full 
?i-point HTL functions have not been consistently given yet Even with such a limited knowledge on the n-point 
HTL functions in the RTF, analyses of the 3-point functions have shown [l^ that there exist the HTL functions 
being proportional to g^T^ (g: the coupling constant, and T: the temperature of the environment), never existing 
in the ITF where all the n-point HTL functions being proportional to g^T^. Here it is worth noting that in the 
RTF the Ward-Takahashi identities have been confirmed to hold only among the ITF counterparts, namely among 
the special linear combinations of n-point HTL functions corresponding to those in the ITF, having not been seen 
anything among those HTL functions not existing in the ITF. 

Recently beginning of the relativistic heavy ion collision experiments at BNL-RHIC has attracted an increasing 
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interest in studying the physics in thermal QCD. Progress of the physics of quark-gluon plasma and the development 
of the calculational framework towards the system in non-(or out-of-) equilibrium have urged us to determine explicitly 
all the HTL's in the RTF with the relations such as the Ward-Takahashi identities they might have satisfied. 

In this paper we determine explicitly in the RTF, especially in the "physical representation" in the closed-time-path 
(CTP), or the Keldysh formalism the HTL contributions to the n-point vertex functions in thermal QCD/QED 

with 71=2, 3, 4 which are in urgent demand. Also studied are the structure of the Ward-Takahashi identities that might 
be satisfied among them. As an application we write down the Dyson-Schwinger equations in the HTL approximation 
for the fermion mass function. Results can be summarized briefly as follows 

i) In the RTF there are two types of n-point vertex functions having the HTL contributions; one is the vertex 
function with n external gauge bosons (ng vertex function) and the other with a pair of external fermions and (n-2) 
external gauge bosons (2f-(n-2)g vertex function). There arc no other n-point HTL functions, thus there are no 
n-point HTL functions with external ghost lines. This situation is completely the same as in the ITF [Q]. 

ii) All the n-point HTL vertex functions with one retarded index in the physical representation in the CTP formalism 
take exactly the same expressions (thus are proportional to g^T"^) as those analytically continued from the n-point 
HTL vertex functions determined in the ITF through simple but corresponding continuation paths. They satisfy 
the same Ward-Takahashi identities as those in the ITF 0, . The n-point vertex functions with one retarded index 
(Vraa, Vara, Vraaa, Varaa, etc) are nothing but the functions exactly corresponding to the n-point vertex functions 
that can be calculated in the ITF, and are given by the specific sum (for details, see section II) of those in the "single 
time representation" in the CTP formalism, where a given n-point Feynman diagram contains a Keldish index at each 
end (external vertex) taking values {1,2}, corresponding to the two branches of the closed-time-path contour. 

iii) The n-point (n=2, 3, 4) vertex functions with two retarded indices (Vraaj Vrraa, etc) are characteristic of the 
RTF, never appearing in the ITF. The corresponding n-point HTL vertex functions must be classified into two groups: 
the ng vertex functions and the 2f-(n-2)g vertex functions. Other types of n-point vertex functions in general never 
have the HTL contributions as mentioned in i) above. It is quite remarkable that in QCD the ng vertex functions 
with two retarded indices are proportional to g^T^ in contrast to the 0{g^T^) behavior of the usual HTL vertex 
functions, thus are expected to play important roles in studying the temperature effects. Nevertheless we can verify 
that they satisfy the QED-type Ward-Takahashi identities between the corresponding ng- and (n-l)g HTL vertex 
functions, still guaranteeing the gauge invariance of the HTL approximation. Thus we can prove, through the explicit 
calculations of the n-point HTL vertex functions, the gauge invariance of the real time thermal QCD/QED in the 
HTL approximation. It is also worth noting that the HTL contributions to the 2f-(n-2)g vertex functions with two 
retarded indices totally vanish , thus in QED there are no additional HTL fimctions other than those appeared in 
the ITF. 

iv) In performing the calculation of the Feynman diagrams in the RTF we should be very careful for the treatment 
of the singular functions (the Dirac-(5 function and the principal part) appearing in the free propagators. We must use 
the properly regularized forms during the calculation and should take the limit e — *■ at the end of all calculations, 
as noted by Landsman and van Wheert Q| . 

v) We explicitly write down the HTL resummed Dyson-Schwinger equation for the physical fermion mass function, 
namely the retarded fermion self-energy function, in thermal QCD/QED, which can be used to investigate the nature 
of the chiral phase transition at finite temperature. Some comment on the double counting problem is also given. 

This paper is organized as follows. In the next section II we give a brief review of the "physical representation" 



in the Keldysh or the CPT formalism. In section III the n-point HTL vertex functions (n=2, 3, 4) are explicitly 
determined in the physical representation. The necessity of the use of regularized form of the singular functions is 
demonstrated. The HTL Ward-Takahashi identities between the four- and three-point functions together with the 
three- and two-point functions are explicitly verified. As an application the Dyson-Schwinger equations in the HTL 
approximation for the physical fermion mass function E/j is derived in section |^. Conclusions and some discussion 



are given in the last section VI , 



II. REAL TIME CLOSED-TIME-PATH OR THE KELDYSH FORMALISM 

We use the closed-time-path (CTP) or the Keldysh formalism ^, |^ of the RTF throughout this paper. In this 
formalism there are two familiar representations, namely by following the terminology of Ref. the "single time" 
representation and the "physical", or the "retarded- advanced" representation. In the single time representation the 
single-particle propagator for free bosons has the 2x2 matrix form 

M)- ("^"(^^ i?12(i^)A 

[d2i{K) D22iK) ) ' 
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^ij (^)i = 1> 2) are given in momentum space by 

D^K) = - . , (II.2a) 

^-(^) = - ^^2^^^ ' (n-2d) 

A ^ — + z£ A ^ — — 

where fsT = (fco,fc), ^ denotes the step fmietion, and the equiUbrium distribution function is given by ^^(fco) = 
l/[exp(|fco|/r) — 1]. For fermions the bare propagator can be also written as the 2x2 matrix form 

«w=(S;(2)IS))- ("■^' 

and 

= {Ji-m)D^,{K) , (II.4a) 

D,j = D,j{nB{ko) ^ -npiko)} , (^,J = 1,2), (II.4b) 

where the Fermi-Dirac distribution is given by npiko) = l/[exp(|A;o|/T') + 1]- Components of these propagators are 
not independent of each other. There is an algebraic identity 

Gil + G22 = G12 + G21, (II-5) 

where G stands for D or S, respectively. 
By an orthogonal transformation 

G = Q-'^GQ, G^QGQ-^; G = DotS , (II.6) 

Q - ^(l-'^z) , (II.7) 
we arrive at the propagator in the physical, or the retarded-advanced representation, 

^^^)-[ DnA{K) DnniK) ) ' ^^^'^^ 

and Doif} (a, (3 = A,R) are 

Daa{K) = 0, (II.9a) 

Dar{K) ( = Da{K)) = -2 -— , (II.9b) 

A"' — m'' — isgii{ko)e 

DnAiK) ( ^ Dn(K)) = — j— , (IL9c) 

A^ — m"^ + zsgn(Ko)e 

DimiK) ( = Dc{K)) = (1 + 2nB{ko))sgn{ko)[DnA{K) - Dar{K)] , (II.Qd) 



for bosons and 



Saa{K) = 0, (Il.lla) 

SAuiK) ( . SAiK)) = - , (11.11b) 

SuAiK) ( . Sum ^ iJ< - , (11.110) 

Srr{K) ( = Sc{K)) = (1 - 2nF{ko))sgn{ko)[SRA{K) - Sar{K)] , (Il.Ud) 
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for fermions, where the last equations in Eqs. (II. 9) and (H-ll) are consequences of the fluctuation-dissipation theorem. 

Propagator (connected 2-point Green function) G and the 1-particle irreducible 2-point vertex function E satisfies 
the Dyson equation {G ~ D or S, and correspondingly E denotes the bosonic or fermionic 2-point vertex, i.e., 
self-energy function) 



EctsG — (73, 

EcTiG = (Ji, 



GctsE = (73, 
GctiE = (Ti. 



(IL12a) 
(II.12b) 



In the above, the relations (II. 1), ( pl.3| ), ( p^I.5| )-(II.8), (11.10) and (11.12) also hold for full propagators and full 2-point 
vertex functions. 



It is also worth noticing here that in the above equations for the free propagators, Eqs. (IL1)-(II.4) and (11.8)- 
(11.11), we do not use the expressions that contain explicitly the singular functions themselves, but rather use the 
functions with the regularization parameter e, e.g.. 



[K^ -m? + is] 



= PP 



1 



[K'^ - m?] 



(11.13) 



where 



PP- = 

X 



'k5{x) 



(II. 14a) 
(II. 14b) 



Keeping e finite (i.e., £ ^ 0) till the end of all calculations then taking the limit e ^ at the end is extremely 
important as we shall explicitly see in the next section. 



III. THE n-POINT (n=2, 3, 4) VERTEX FUNCTIONS IN QCD/QED IN THE HTL APPROXIMATION 

In this section we calculate the HTL contributions to the one particle irreducible n-point vertex functions (n=2, 3, 
4) in the physical or the retarded-advanced representation. One can construct in general the n-point vertex functions 
in the physical or the retarded-advanced representation from the components of the real-time n-point functions in 
the single time representation as in literature ^, ||, |^, |l3|, Q , in which some care must be taken because the notation 
varies from literature to literature. Throughout this paper we use the notation as follows; a given one particle 
irreducible n-point Feynman diagram in the single time representation contains a Keldish index at each end (external 
vertex) taking values {1, 2}, corresponding to the two branches of the closed-time-path contour, while in the physical 
representation it contains a retarded-advanced index at each end taking values {i?, A], corresponding to the retarded 
or the advanced prescription. Our notation follows the one in Ref. 0. 



A. The n-point vertex functions (n=2, 3, 4) in the physical representation 



For convenience we here present explicitly the 2-, 3-, and 4-point vertex functions in the physical representation 
constructed from the components in the single time representation. Generally speaking the n-point function has 2" 
components. These components obey one constraint equation, which reduces the number of independent components 
to 2" — 1. In the physical representation this is express ed by t he fact that t he n-p oint vertex function with n- 
external advanced indices automatically vanishes, Eqs. (III. la), (111.3a) and (III.4a), below. In equilibrium, the 
Kubo-Martin-Schwinger conditions impose additional constraints, reducing the number of independent components 
to 2^"^^^ — 1. For details of the construction, see Refs.Q, Since the transformation formula between the vertex 
functions in two representations do not take care of the external particle species (fermion or gauge boson) , here we 
simply write Vai3...5, {a, /?, S = R, A}, or Vij,,,i, {i,j, ■■■,1 — 1, 2} for vertex functions in the physical and the single 
time representations, respectively. 
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1. 2-point vertex functions 



Vaa = 0, (Ill.la) 

Vra {^Vr) - Vii + Vi2, (Ill.lb) 

Var {^Va) = Fii+l/21, (III.lc) 

Vrr ( = Vb) = 1^1+1^22. (Ill.ld) 

Among them the following relations hold, 

Var{K) = Vra*{K), (III.2a) 

Vrr{K) = {l + T^2n{ka)}sgxi{ko){VRA(K)^VAR{K)}, (III.2b) 

where rj = 1) for boson (fermion) and n(fco) is the corresponding equilibrium distribution function. 

2. 3-point vertex functions 



Vaaa 


= 0, 








(111.3a) 


Vraa 


= ■^111- 


~ V112 - 


1- ^121 - 


- V122, 


(III.3b) 


Vara 


= Viii - 


' V112 ' 


1- V2n - 


^V212, 


(IIL3c) 


Vaar 


= Viii- 


- V121 - 


- V211 - 


- V22I, 


(III.3d) 


Vrra 


= Viii- 


- V112 - 


" ^221 " 


" V222 , 


(III.3e) 


Vrar 


= Viii- 


- V121 - 


' ^212 " 


- V222 , 


(III.3f) 


Varr 


= Viii- 


~ V211 - 


" V\22 " 


- V222 , 


(IIL3g) 


Vrrr 


— Viii - 


~ ^122 " 


V V212 - 


1- V^221. 


(III.3h) 



The Kubo-Martin-Schwinger conditions impose additional 4 constraints, reducing the number of independent compo- 
nents to 3. 



3. 4^-point vertex functions 



Vaaaa 
Vraaa 
Varaa 
Vaara 
Vaaar 

Vrraa 

Vrara 
Vraar 



0, 

1^111 
1^111 
1^111 
1^111 
1^1111 

Villi 



V1112 
V1112 
V1112 

^1121 

V1112 
V1112 
V1121 



^1121 
^1121 
^1211 
^1211 
^1121 



^1211 

V2111 
V2111 
V2111 

^1122 

V12II + V1212 + V212I 
V12II + V122I + V2112 



^1122 
V1122 
V1212 
Vl22\ 
V22II 



^1212 - 


- V122I - 


- V1222 


V2II2 - 


- V212I - 


- V2122 


V2II2 - 


- V22II - 


- V22I2 


V212I - 


- V22II - 


^ ^2221 


V22I2 - 


- V222I - 


- ^2222 


V2122 - 


- V222I - 


- V2222 


V2122 - 


- V22I2 - 


- V2222 



(III.4a) 
(IIL4b) 
(III.4c) 
(III.4d) 
(III.4e) 
(III.4f) 
(III.4g) 
(IIL4h) 



There are seven independent vertex functions, Eqs. (III.4a-h), in this case, and the other eight vertex functions can 
be obtained from them using the Kubo-Martin-Schwinger conditions. 



B. The n-point vertex functions (n=2, 3, 4) in the HTL approximation 

Now we calculate the n-point vertex functions (n=2, 3, 4) in the physical representation in the HTL approximation. 
In the following we consider the massless QCD/QED, namely the high temperature hot QCD/QED, where all the 
fermions and the gauge bosons are massless. 
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Because we are accustomed to the Feynman rules in the single time representation, we calculate explicitly the 
r ight- hand-side of Eqs. (III.l), (III. 3) and (III.4). As we have stressed at the end of the last section we use Eqs. 
([1.13) and (11.14) for propagators, namely, we do not use the expressions that contain explicitly the singular functions 
themselves, but rather use the functions with the regularization parameter e and keep e finite (i.e., e ^ 0) till the end 
of all calculations then take the limit e — > at the end. 



1. 2-point vertex functions, or the fermion self-energy and the gauge boson polarization tensor 

The 2-point vertex function in QCD/QED is usually named as the self-energy part E for fermion, while it is called 
as the (vacuum) polarization tensor H^" for gauge boson (gluon/photon). Although the HTL results of E and H^'^ in 
the single time representation as well as in the ITF are already well-known ||^, |l^, |l^ , for the sake of completeness we 
here present the results in the physical representation. In fact in calculating their HTL contributions in the physical 
representation, some care being worth mentioning must be taken, obtaining the result that could not be seen, at least 
in the explicit form, in the previous calculations ^ p^ . 

The HTL contribution to the fermion self energy, ST,, in the single time representation is obtained in QCD by 
calculating the diagram shown in Fig.|l| (in QED g^Cp should read e^); 

^D,^{K + Q)D,,{K), P + Q^O, (111.5) 

where and Dij are given in Eqs. (II. 2) and (II. 4). With the use of Eqs. (III.l) we get the HTL fermion self energy 
in the physical representation. There arc nothing new in the obtained results in QCD, thus only reproducing the 
previous results; 



di:RA{P,Q) - -^T' / dn—4- , (in.6a) 

o^TT J Q ■ K + ie 
SJ^ARiP^Q) - ~^^T^ I dil ^— , (IIL6b) 



327r J Q-K-ie 
5^rr{P,Q) - 0, (IIL6c) 



where = (1, k) with k = k/k, k ~ vfc^, being the unit three vector along the direction of k. 

As was expressed in Eqs. (III.l) we usually denote the {i?A}-component as the R-, or the retarded component, and 
the {Ai?}-component as the A-, or the advanced component of the corresponding quantity. Namely E/j^ = S/j is 
the fermion self-energy part of the inverse retarded fermion propagator, having a definite physical meaning, i.e., the 
physical fermion mass function. 

As for the HTL contribution to the gauge boson polarization tensor, SH^'^ , the diagrams to be calculated in the 
single time representation are shown in Fig. || (P + Q = 0) ; 



-4e2(-l)«+^-2 / AA^ [2if^if- - g^^'^K^] D,,{K)Dj,{K + P), (QED) 

.14 (HI. 7) 

252(_l)i+J-2 J [2K^K'' - gt^'^K'^] 

X {NcD,,{K)Dj,{K + P) - NfD^j{K)Dj^{K + P)} , (QCD) 



from which we get the HTL gluon polarization tensor in QCD (for the photon polarization tensor in QED, g'^{Nc 
{l/2)Nf) should be understood as e^). 



(P, Q) = (a^c + In,) J dn {K^K^^f^^ - 5^°.^°) , (ni.8a) 
(p, Q)--'-^{n. + 1a.,) / dn {K^K^jf— - s^V') , (ni.8b) 

S^^niP^ ^) - (^^ + l^f) J - ^) • (ni-8c) 

Same as the fermion self-energy function ^jia^^ab) = ^li (^'a) gauge boson polarization tensor in the inverse 

retarded (advanced) gauge boson propagator, thus having a definite physical meaning. 
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It should be noted that the {i?i?}-component of the polarization tensor, SIl'^j^, is proportional to g^T^, compared 
to the ordinary retarded or advanced polarization tensor, '511^'^ or '511^'^, being proportional to g^T^ ■ Also noted 
is that all the HTL contribution in the ITF are of the order 0{g^T^), thus the above results show completely new 
vertex functions appear in the RTF. 



2. 3-point vertex functions 

There are two types of 3-point vertex functions, the fermion-gauge-boson vertex function and the 3-gauge-boson 
vertex function. Needless to say in QED only the fermion-gauge-boson vertex function exists. 

a. The fermion-gauge-boson vertex functions We define the fermion-gauge-boson vertex function in the HTL 
approximation as in Figj^, where 

*r'^ = jf + ST", (III.9a) 

with 7^ representing the tree vertex 

= , otherwise, (III. 9b) 

and the HTL contribution to the fermion-gauge-boson vertex function, JF^, is obtained by calculating the one-loop 
diagrams shown in Fig.^ 

\ ^^e^ J Ak^' , (QED) 

^rf,jP,g,i?) = <^ ^ ~ ^irr ,~ \ (IILlOa) 

I ^ig'Cp J - 2ig^N, J K [v^k + V^k) , (QCD) 

where, 

V,,k = {-lT+'+''-^D,j{K)Djk{K-Q)Dk^iK + P) , (IILlOb) 
V^Jk = {~ly+'+''-^D^,{K)DJk{K~Q)Dk^{K + P) . (IILlOc) 

The results in QCD in the physical representation are given as follows (in QED g'^Cp should be understood as e^); 

o^TT J {K-P - ie){K-Q + le) 

a'^T^ f kf^ k 

ST%AP,Q,R) = -%t-Cf / dn 2 , (Hl.llb) 

oiTT J (K-P -\- ie){K-Q - le) 

a'^T'^ f k^ k 

•srSsjWO.fl) = iri;js(P,o,ii) = ir^j^iP.Q.R) = «r"„^(p,(3,R) = o . (iiLiie) 

We can see that all the HTL terms of the fermion-gauge-boson vertex function, SV^ , are proportional to g^T^ . 

In obtaining above results it is importa nt t o use the re gulariz ed expressions for free propagators with the finite 
regularization parameter e ( ^ 0), Eqs. (II.1)-(II.4) and (II.13)-(II.14), not using the expressions in terms of the 
explicit singular func tions fr om the beginning. To see this point more clearly l et us calculate the above HTL vertex 
function (5r^^^, Eq.([II.lla), explicitly in QED. In each calculation of (Ill.lla-d) we must evaluate four diagrams in 
the single time representation, Eqs. (III. 3), thus face the calculation, e.g., oi 5V'^^^{P,Q,F(), 

SV^^P, Q, R) = 4ie^ J ^D,,{K)DniK - Q)D,,{K + P). (IH.12) 
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The loop- momentum integration over if should be performed by keeping every Si finite (i.e., Si ^ 0). Then the 
singularities of the integrand as a function of fcp only poles in the complex /cQ-plane thus the integration over 
ko can be carried out through the residue analysis. By neglecting the 0{e^T) contributions we get after some 
manipulations (here we set £i = £2 = £3 = e) 



5r5'ii(p,Q,i?) = - 

where n(fc) — nsik) + np{k) and = 



(27r)3 



(27r)3 

„2 



27r)3 



dkkn{k)dnKl^ 



dkkn{k)dnK^ ^ 



dkkn{k)dnKl^^ 



K+Q 





h (£) 


- k/k. Similar calculation [ 


jives 


£ e 




+ {K+PY + {sY {K+QY^ 


HsY 


K+P e 




+ {K+Pf + [sY [K+QY 


+ ieY 


£ K+Q 




'{K+PY + {£Y {K+QY + 


(sY ■ 



(III. 13) 

(III.14a) 
(111.14b) 
(IIL14c) 



Adding Eqs.( |III.13D , ( [III. 144 c) we get Eq.( |[II.llaD . 

If we take naively the limit ^ before the loop-momentum integration over K^^, and use in Eg. ( [II. 12 ) the free 
propagators explicitl y cont aining the singular Diiac-S function and the principal part, then the naive manipulation 
gives, instead of Eg. ( [II. 13 ), 



ST'^,,{P,Q,R) 



(27r) 



dkkn{k)dnK^'^^ 



PP- 



1 



tPP- 



1 



K+P K+Q 



-5{K+P)5{K+Q) 



(111.13') 



Other three vertex functions coincide with the limit of Egs.( III.14a -c) and we can not get Eg. ( III. 11a ). 

Additional term in Eg. (III. 13') being proportional to the product of Dirac-(5 functions may have its origin from the 
integral whose integrand being the product of two principal parts 



dy / dxPP[- 



1 



a{y) 



IPPf 



1 



KvY 



(III.15) 



If, at some y = yo inside the integration range over y, a{y) = b(y) happens then we face the integration over x of the 
sguare of principal part, which can not be well-defined. This is exactly what happens in the above calculation of Eg. 
(111.12) with the free propagators explicitly containing the singular Dirac-(5 function and the principal part. Without 
concrete and consistent prescriptions how to treat the product of singular functions we can not get a definite result, 
and the e-regularization method gives such a prescription. As we have already noted throughout this paper we use 
the £-regularized singular functions. 



b. Three gauge-boson (three gluon) vertex functions We define the three gluon vertex function in the HTL ap- 
proximation as in Fig.^, where 

*Yt,up ^ yt^up ^pM-^P^ (III.16) 

with V^^"^^ representing the tree vertex, which does not exist when the number of external retarded indices are zero 
or two. The HTL contribution to the three gluon vertex function, JP^"'', is obtained by calculating the one- loop 
diagrams shown in Fig.^ 

5rf;/(P,Q,i?) = - '-^{-ly+^+i-^ j d^KKt^K-KP 

X [N,D,j{K)Djk[K - Q)Dk^{K + P) - Nf^,^{K)D,t{K - Q)Dk^{K + P)] . (III. 17) 
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By using Eqs.(III.3) we get the HTL results in the physical representation; 



9^ 
127r 



1 



127r 



5r7^''^(P,Q,i?) = - 



5^ 
127r 



5^ 
127r 



67r 



67r 



+ 2^/ 



dflK^'K^RP 



dnX'^K'^KP 



dnXf^K^KP 



{K-P - ie){K-R + ie) [K-Q + ie)(K-R + ie) 

(111.18a) 

90 

{k-P + ie){k-R + ie) (i^T-g - ie){k-R + ie) 

(111.18b) 

Po go 

(iiT-P + ie){k-R - ie) {k-Q + ie){k-R - ie) 

(III.lSc) 

Po 90 

{k-P - ie){k-R - ie) {k-Q - ie){k-R - ie) 

(III.18d) 

1 1 

{k-P - ie){k-Q + ie) {k-P + ie){k-Q - ie) 

(Ill.lSe) 

1 1 

{k-P + ie){k-R - ie) {k-P - ie){k-R + ie) 

(III.18f) 

1 1 



{K-Q - ie){K-R + ie) {K-Q + ie){K-R - ie) 

(Ill.lSg) 
(III.18h) 



It should be noted that depending on the number of external retarded indices the HTL three gluon vertex functions 
can be classified into two groups; i) number of i? is 1 or 3, and ii) number of i? is or 2. Any vertex belonging to the 
first group has the tree vertex and is proportional to g^T^, while other vertices belonging to the second group do not 
have the tree terms and are proportional to g^T^, having no counterparts in the ITF. 

3. 4-point vertex functions 

There are two types of 4-point vertex functions, the fermion pair and 2-gauge boson vertex functions and the 
4-gauge boson vertex functions. In QED, as wc shall see below, all the HTL contributions to the 4-photon vertex 
functions vanish. 

a. One fermion pair and 2-gauge boson vertex functions We define the fermion pair-2-gauge boson vertex func- 
tions in the HTL approximation as in Fig.^, where 



(IH.19) 
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and the HTL contribution to the fermion pair-2-gauge boson vertex function, ^r^*^, is obtained by calculating the 
one-loop diagrams shown in Fig.^ and their exchanged diagrams between the external legs with momenta R and S\ 

-i8e2(_l)»+J+fe+' / ^Kr^K'^ ^D,,iK)D,k{K - Q)Dm{K + P + S)Du{K + P), 

(QED) 

^i8.g2CF(-l)'+^+'^+' / J^^K^K" ^ 

X [CFp,,iK)D,kiK - Q)Dki{K + P + S)Du{K + P) 
-N,D,,{K)D,k(K - Q)Dki{K + P + S)Dh{K + P) 
+ \n,D,,{K + P)Dji{K - P - S)Dik{K + R)Dk^{K)\^ ■ (QCD) 

(III.20) 

The results in QCD in the physical representation are, with the use of Eqs.(III.4), given as follows (in QED g^Cp 
should be understood as e^); 



st'X'aar(p,Q^R,s) - 



S^RRAAiP^Q^R^S) 



9^ 
327r 



Cf / dnRf^K" fi 



1 



1 



K-P -ie K-Q + ie 



1 



+ 



1 



KiP + S)-ie K-{P + R)-ie 



(III.21a) 



9^ 
327r 



Cf / dnkf^K" ^ 



1 



1 



K P + ie K-Q - is 



K-{P + S)+ie K-{P + R) + ie 



(III.21b) 



327r 



9^ 
327r 



Cf / dnK^'K" fi 



1 



K-P + ie K-Q + ie 



1 



1 



K-{P + S)+ie K-(P + R)-ie 



(IIL21c) 



Cf / dnRf^K" fi 



1 



1 



K P + ie K-Q + ie 



1 



+ 



1 



K-{P + S)-ie K-{P + R) + ie 



^^^ARAiP^ R^ s) = Sr^RAARiP^ Q, R, S) - 0. 



(III.21d) 
(IIL21e) 



Other eight vertex functions can be obtained from Eqs. ([II.21a-e) using the KMS conditions. Again all the HTL 
contributions to the verte x funct ions with external fermion legs are of 0{g^T^). 

Some comment on Eg. ( III. 21^ ) must be given. Hou Defu et al, Ref.||T^, calculated the same HTL vertex functions 
and claimed the existence of the non-zero HTL contributions, only the first function ^r^'^^^(P, Q, R, S) to vanish in 
equilibrium. However, as we explicitly show, all the three functions vanish in the HTL approximation so long as the 
propagation of thermal " quasi-particles" can be described by the form of free particle propagators, Eqs. ([L1)-(II.4), 
no matter which in equilibrium nor just out-of-equilibrium. 



b. ^-gauge boson vertex functions The HTL contribution to the 4-photon vertex function in QED completely 
vanishes, thus we confine our interest to the 4-gluon vertex functions in QCD. Because of the complexity in their 
tensor structure and of our interest in their application, we calculate the HTL's for the vertex functions being summed 
in the color indices over two of the external gluon legs. 

Defining the 4-gluon vertex functions in the HTL approximation as in Fig||, where the summation over the color 
indices of the gluons having momenta R and S should be understood and 

*Yi^i'P<y ^ ^^t^i'P'y ST'^'P", (HL22a) 

with M/^"'"^ representing the tree vertex, 

w^X = (-r'(W'^-3'^v^-.9'^v), ^=j = k = i, 

= 0, otherwise. (HL22b) 

The HTL contribution to the 4-gluon vertex function, ST'^'^''"' , is obtained by calculating the one- loop diagrams shown 
in Fig.ll^ and their all possible exchanged diagrams among the external legs with momenta Q, R and S (Note that 
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the color indices of the gluons with momenta R and S are summed over) ; 



, 32 



[Nl {V,,ki{P. Q, R, S) + Kj7fc(P, Q, S, R)} - CfNf {v^JkliP, Q, R, S) + V,jik{P, Q, S, i?)} 



+^V,kji {P, R, Q, S) + ^V^kji (P, R, Q, S) 



where, 



V.,ki = {-iy+'+''+'D,j{K)D,k{K - Q)DkiiK + P + S)Dh{K + P), 
V.jki = {-iy+'+''+'D,j{K)Djk{K-Q)DM{K + P + S)Du{K + P). 
The resuhs in QCD in the physical representation are, with the use of Eqs.(III.4), given as follows; 



S^'nTAAiP^Q^R^S) 



ST'i,JnAiP^Q,R,S) 



9^ 



Po + ro 



dnK>'K''KPK^ 

go 



1 



1 



K-S + ie K-R + ie 



K-{P + R)-ie K-Q-ie K-P~ie K-{Q + R)+ie 



Po + ro 



1 



1 



K-S + ie K-R + ie 
Po , qo+ rn 



sr'X7RAiP^Q,R,s) = 



K-{P + R) + ie K-Q^is K-P + ie K-{Q + R)-ie 



6tt 



Po + ro 



dVLKPK^KPK^ 

qo 



1 



1 



K-S + ie K-R-ie 
Po , qo + ro 



K-{P + R)-ie K-Q + ie K-P + ie K-{Q + R)-ie 



g^T^ 



6tt 



Po + ro 



dnKf^K^K'^K" 

qo 



1 



1 



K-S-ieK-R + ie 
Po , qo+ ro 



K-{P + R)+ie K-Q + ie K-P + ie K-{Q + R) + ie 



2rp3 



67r 



Nc + \Nf 



dnK''K''KPK^{P ~Q)-K 
1 1 1 



K-{P + R)+ ie K-P + ie K-{Q + R) ~ ie K-Q - ie 



c.c 



2rp3 



6tt 



1 



1 



dVlK'^K^K^K" 
1 



K-P + ie K-{Q + R)-ie K-R - ie 



^KTAB.iP^Q,R,s) = -■- 



2rp3. 



Nc + ^Nf 



dnK'^K^K^K'^ 
1 



K-P - ie K-{P + R)-ie K-S + ie 



c.c 



(III.23a) 



(III.23b) 
(IIL23c) 



(III.24a) 



(III.24b) 



(III.24c) 



(III.24d) 



(IIL24e) 



(III.24f) 



(IIL24g) 



Other eight vertex functions can be obtained from Eqs. (III.24a-g) using the KMS conditions. It is worth noting that 
in this case again depending on the number of external retarded indices the HTL four gluon vertex functions can be 
classified into two groups; i) number of i? is 1 or 3, and ii) number of i? is or 2. Any vertex belonging to the first 
group is of 0{g^T^), while other vertices belonging to the second group are are of 0{g'^T^), having no counterparts 
in the ITF. The existence of 0{g^T^) HTL contribution is characteristic of the n-gauge boson vertex functions. 
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IV. WARD-TAKAHASHI IDENTITIES IN THE HTL APPROXIMATION 



Having determined the 2-, 3-, and 4-point vertex functions in QCD and in QED in the HTL approximation, in this 
section we verify that the Ward-Takahashi identities satisfied between them are all of the simple QED-type identities. 
As we have shown the HTL n-gauge boson vertex functions with even number external retarded indices show, con- 
trasting to the ordinary 0{g^T^) behavior of other vertex functions, the 0{g^T^) behavior being totally absent in any 
amplitude in the ITF. We here verify that they also satisfy the simple QED-type Ward-Takahashi identities between 
themselves in the HTL approximation, thus verifying with an explicit calculation the gauge invariance of thermal 
QCD/QED in the HTL approximation in the RTF, being guaranteed by the QED-type Ward-Takahashi identities. 



A. Ward-Takahashi identities between the three- and four-point vertex functions in QCD/QED 

1. Between the 3- and 4-point HTL vertex functions with a fermion pair legs 

There are no differences between the 3- and 4-point HTL vertex functions with a fermion pair legs in QCD and in 
QED except the group factor, thus we only study in QCD. Comparing the results Eqs.(HL21) and Ec|s.(III.ll) we 
obtain the following Ward-Takahashi identities between vertices shown in Fig.^ (i.e., ||) and FigJ| (i.e., 

*^RAAAiP^ s) = ST'^j.^AP + R, Q, s) - sv^n^AP, Q + R, s), (iv.ia) 
Rt^ *^ARAAiP^ R^ s) = ^^ARAiP + R. s) ~ 5v\^a{p^ Q + R, s). (iv.ib) 

R^^ *'^aara{P^ R^ S) = 5T'i,AA{P + R. Q, S) - ST'^^j.^P, Q + R, S), (IV.lc) 
*^AAARiP^ R' S) = SV^^jiiP + R, Q, S) - ST'XARiP, Q + R, S). (IV.ld) 
When we operate in place of i?^, then obviously the role of S and R are exchanged in Eqs.(IV.l) and also the role 



of Eqs.( IV.lc) and (IV.ld) should be exchanged. 



2. Between the 



irf 4-gluon HTL vertex functions 



In QED there are no 3-photon vertex functions and the 4-photon vertex functions totally vanish in the HTL 
approximation, thus we study in QCD the relation between the 3- and 4-gluon HTL vertex functions. Comparing 
the results Ens. (III. 24) and Ec|s_.(III.18) we obtain the following Ward-Takahashi identities between vertices shown in 
Fig.| (i.e., |l^) and Fig.| (i.e. 



^fJ."'' RAAA\^' ^ 


^,R,S) 


= S^'AAiP^ 


^R,Q,S)- 


S^R'AAiPyQ- 


hi?, S), 


(IV.2a) 


Rp-^^ARAAiP^ ^ 


^,R,S) 


= ST'^/i^AiP^ 


hi?,g,5)- 


sr7RAiPyQ- 


^R, S), 


(IV. 2b) 


Rp-^^AARAiPy ^ 


^,R,S) 


= sr'g'up^ 


hR,Q,S)- 


sr7RAiPyQ~ 


^R, S), 


(IV.2c) 


Rt^^^AAAR(Py ^ 


^,R,S) 


= sr7lRiP^ 


hR,Q,S)- 


sr7lRiP,Q~ 


^R, S), 


(IV. 2d) 


Rt^^^RRAAiPy ^ 


^,R,S) 


= sn'^AiP^ 


hR,Q,S)- 


sr^RAiP^Q- 


^R,S), 


(IV.2e) 


^fJ."^ RARA\^y ^ 


^,R,S) 


= -srj{';^APyQ + R,s), 






(IV.2f) 


Rf^^^RAARiPy 


^,R,S) 


- sr'^'UP^ 


hR,Q,S)~ 


sr^ARiP^Q- 


\'R, S). 


(IV.2g) 



When we ope rate S'ujii_place of i?^, then obviously the role of S and R are ex change d in E qs.(IV .2) and also the role 
of Eqs.([V.21) and (IV. 2g) should be exchanged together with the role of Eqs.(IV.2c) and (IV. 2d). The situation is a 
bit complicated in this case and the resulting identities are explicitly reproduced: 



^ fJ'^^RAAAiP^ R^ ^) 

S^6r7^^^^{P,Q,R,S) 
S^5T7rRA{P,Q,R,S) 

^H^^AAARiP^ R' ^) 
^ fJ^^^RRAA(P' R' ^) 
^ fJ^^^RARA(P' R' ^) 

s.sr'^RrARiP^Q^R^s) 



sr7lAiP,Q 
sr7RAiP^Q 



(IV.Sa) 
(IV.3b) 
(IV.3c) 
(IV.3d) 
(IV.3e) 
(IV.3f) 
(IV.3g) 

It is worth noting that Eqs.( IV.2^ g) and (IV.3e-g) are the new QED-type Ward-Takahashi identities satisfied among 
the vertices with the 0{g^T'^) behavior thus being absent in the ITF. 



sr-^'UP^ 
sr7RAiP 
^n'lRiP^ 
sr^i^AiP^ 

^^R^RAiP 

ST^k'XRiP 
-ST^^AP^Q 



S,Q,R) 
S,Q,R) 

S,Q,R)-6T7^j,{P,Q 
S,Q,R)-5T7^j,{P,Q 
S,Q,R)-ST7^^{P,Q 

S,Q,R)-ST7Xr{P,Q 
S,R). 



S, R), 
S, R), 
S,R), 
S,R), 
S,R), 
S,R), 



13 



B. Ward-Takahashi identities between the two- and three-point vertex functions in QCD/QED 

1. Between the 2- and 3-point HTL vertex functions with a fermion pair legs 



Since only the difFerence between in QCD and in QED is the group factor, by comparing the results Eqs.(III.ll) 
and Eqs.(III.6) it is easy to see that we can verify the same Ward-Takahashi identities in QCD and in QED between 
vertices shown in Fig.|| (i.e.,||) and Figj^: 



(P,Q,i?) 



S1:ra{P, Q + R)- 5^ra{P + R,Q), 
6j:ar{p, Q + R)- SJiABiP + R,Q), 
5^AB.{P, Q + R)- SJ:ra{P + R,Q), 
6^ra{P, Q + R)- 6J:ab{P + R,Q). 



(IV.4a) 
(IV.4b) 
(IV.4c) 
(IV.4d) 



2. Between the 2- and 3-gluon HTL vertex functions 



In QED there are no 3-photon vertex functions, thus we study in QCD the relation between the 2- and 3-gluon 
HTL vertex functions. Comparing the results Eqs.(III.18) and Eqs.(III.8) we obtain the following Ward-Takahashi 
identities between vertices shown in Fig|| (i.e., |^) and Fig.^: 



^fJ."^ RAA 
^fJ.''^ ARA 

-"■m"^ rra 

^n"^ RAR 

ARR 

-"-m"^ rrr 



{P,Q,R) 
{P,Q,R) 
{P,Q,R) 
{P,Q,R) 
{P,Q,R} 
{P,Q,R) 
{P,Q,R) 



STIrWQ 
^KWQ 
<5n7^(P,Q 
SU'g'^{P,Q 
SK%{P,Q 
-SIlR^RiP^ 

sn^jl^iP^Q 



SK'AiP 

5^7r{p 



-R) 
-R) 
-R) 

-R), 
R,Q), 

-R)-sn''UP 



R,Q), 
R,Q), 
R,Q), 
R,Q), 



R,Q). 



(IV.5a) 
(IV.5b) 
(IV.5c) 
(IV.5d) 
(IV.5e) 
(IV.5f) 
(IV.5g) 



Here again we should note that Eqs.( fv\5e| -f) are the new QED-type Ward-Takahashi identities satisfied among vertex 
functions with the Ojg^T^) behavior thus being absent in any amplitude in the ITF. 

In the last section [II we showed that in the physical representation in the real time CPT formalism there exist those 



vertex functions with two (presumably even number in the arbitrary n-point case) external retarded indices having the 
high temperature behavior of 0{(^T^). In the ITF there are no n-point functions with the 0{g^T^) behavior, namely 
all the HTL's are of 0{g^T^) [Q, among which the QED-type Ward-Takahashi identities are satisfied, guaranteeing 
the gauge invariance of the HTL's. In this section we verify that those vertex functions with the high temperature 
behavior of 0{g^T^) also satisfy among themselves the simple QED-type Ward-Takahashi identities in the HTL 
approximation, thus can show explicitly the gauge invariance of the HTL's in the real time thermal QCD/QED. 



To close this section it is better to make notice on the identities ([V.21) a nd (IV.Sg), and also on those (|IV.5e ) and 
( IV.Sl ). They seem to have a bit different structure from others; identities ( IV.21 )77 [V-3g ), ( IV. 5c ) and ( IV. 51 ) have 
the right-hand-sides with a single HTL vertex functions, while the right-hand-sides of other identities in Eqs.(IV.2), 
(IV. 3) and (IV. 5) consist of difference of two vertex functions, b eing fa m iliar in t he Wa rd-Ta kahash i identities in 
QED. We should note, however, that the right-hand-sides of Eqs.( |IV.2| ), ( |IV.3gD , ( |IV.5c) ) and ( |IV.5fD are the HTL 
contributions to the 3-, or 2-gluon vertex functions with two retarded indices, which are e ssentia lly consist of the 
differen ce of two vertex functions. With this fact we can understand that all the identities ( lV.2a -g), ( IV.3a -g) and 
([V.5a-g) actually have the same structure. Further discussion on this point will be given in the last section VI. 



V. DYSON-SCHWINGER EQUATION IN THE HTL APPROXIMATION 



Having determined all possible ingredients, namely the HTL contributions to the 2-, 3-, and 4-point vertex functions, 
now we can write down the Dyson-Schwinger (DS) equations in the HTL approximation. Before doing this, however, 
let us here remember the reason why we have calculated the HTL vertex functions in the physical representation in 
the real time CTP formalism of thermal field theory. 

As we have already noticed in section 111, to investigate the consequences on the physical mass we need to study 



the mass function of the inverse of the retarded propagator. To investigate the chiral phase transition we need the DS 



14 



equation for the retarded component of the fermion mass function, E/j, and to study the magnetic screening the DS 
equation for the retarded component of the gluon polarization tensor, H'^ ■ Nevertheless in many analyses the 
{ll}-component of the fermion self energy in the single-time representation, Sn, not the Eij, have been studied by 
neglecting its imaginary part without much attention. As is well known, ReTij^ = i^eEn, but not for their imaginary 
parts. The lesson from the HTL resummcd perturbation theory jl^ , remainds us of the fact that the imaginary 
parts of Efl and of 11^'' really get the important thermal effects, which in turn affecting their real parts. In this sense 
it is important to correctly construct the DS equations exactly for S^j and for ■ 

Now we are ready to write down the DS equations. For definiteness in this paper we only give the DS equation for 
the fermion self-energy in the HTL approximation, that can be obtained by applying the following approximation 
to the full DS equation; 

i) replace the full gauge boson propagator with the HTL resummed propagator, and 

ii) approximate the full vertex functions to the HTL resummed vertex functions. 

Then we get in QCD the desired DS equation (in case of QED, g^Cp in Eq.(VJ) should read and in Eqs.(V.3) 
should be properly understood), 

- = — y w 

X mAAi-P^ K,P- K)Sra{K) *T''nAAi-K, P,K-P) 
2T 



X ■ 



-[*GraMP-K)- *GarMP-K)] 



Pa - ko 

+ *^RAAi'P,K,P-K)(l - 2nF{ko))sgn{ko)[SRA{K) - SABiK)] 

X *nARi~K, P,K~P) *GraMP - K)} , (V.l) 

where *r'' is the HTL resummed quark-gluon vertex, defined in Eqs.(III.9) and given in Eqs.(III.lO) with Jr^a = 
^AAR ~ 7^- is the HTL resummed gluon propagator, where retarded/advanced propagator is given by 

*G't\,,J-K,K) = ;t7i A^"' + B^"' 

" ' " *H^/^(i^)-i^2^isgn(fco)e *Tll' ^ [K) - T i^g^{ko)e 

^ (V.2a) 



K"^ =F *sgn(fco)£ 

where is the gauge fixing-parameter = in the Landau gauge) and the {_Ri?}-component by 

IT 

*G>i^^{-K,K) = — { *G%^^{-K,K) - *G^/^{-K,K)} , (V.2b) 

Kq 

with *H;^^'* and *H^^'* being the HTL contributions to the retarded/advanced gluon self-energy of the transverse 
and longitudinal modes, respectively. 



.nr(AO^-»?M_i!(i + |io.^), (V,3., 



2^ 



where nig denotes the thermal gluon mass (in QED, g'^{Nc + Nf/2) should read e 

In Eq.( |V.2£j ), A^"',B^"' and D^"' are the projection tensors |jl|| 

A^"'{K) EE gf"" - B'"'{K)- D'"'{K), (V.4a) 
B'^iK) -- (V.4b) 

D^%K) ^ (V.4c) 

k^' EE (fc, fcofc), = -K^ = -{kl - fc2), (V.4d) 
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where k = k/k is the unit three vector along the direction of k. 
Sra/ar is the retarded/advanced fuU fermion propagator, 



Sra/ar{—P, P) — 



and Srr, the {_Ri?}-component or the correlation, 

Srr{-P,P) = (1 - 2r^^.(po))sgn(po){5i^A(-P,P) - Sar{-P.P)}. 



(V.5a) 



(V.5b) 



In the present HTL approximation, the DS equation (V.l) becomes an integral equa tion for the unknown fermion 
self-energy function Y^r. It is worth giving some comments on this DS equation, (V.l): 

i) The HTL resummcd quark-giuon vertex function is substituted for both of vertices in Eq.(V.l). This may cause 
the double counting problem when the loop- momentum becomes " hard" . Thus in the actual analysis we should 
introduce an intermediate momentum scale to cut the loop-momentum integration. In the hard loop- m omen tum 




11| are totally and exactly accounted 
Inclusion of the contribution from 



region the naive ladder approximation may work, while in the soft loop-momentum region we must use Eq.(V.l) 

ii) There are no contributions from the fermion-pair-2-gluon vertex function in Eq.(V.l). In QCD this happens quite 
miraculously, namely the contribution to Y,r from the diagrams shown in Fi 
for through the loop diagram with two 3-point quark-gluon vertices in Eq.(| 
Fig.|ri| causes the trouble of double counting of diagrams. 

iii) The first term in the curly bracket, which exists also in the limit of ladder approximation, has been dismissed 
in previous DS equation analyses | p6[ . As can be seen this term could produce a n do minant contribution due to the 
presence of the enhanced temperature dependence. Result of the analysis of Eq.( V.l) will be given elsewhere [20[|. 



Similarly we can write down the DS equation for the gauge boson polarization tensor, with which the full QCD 
analysis may be performed. In the DS equation for the gluon polarization tensor, which is not given here explicitly, 
the contribution from the gluon 4-point vertex does exist. However, in this case also contributions from the HTL's, 
jp/ij/pcr^ Eqs.(III.24), are totally and exactly accounted for through the loop diagram with two 3-point gluon vertices, 
thus to avoid the double counting problem only the contribution from the tree 3-gluon vertex survives. Results of 
the analysis including the DS equation for the gauge boson polarization tensor itself will also be presented elsewhere 



VI. CONCLUSIONS AND DISCUSSION 



In this paper we calculated the n-point HTL vertex functions in QCD/QED for n= 2, 3 and 4 in the physical 
representation in the RTF. The result showed that the n-point HTL vertex functions can be classified into two 
groups, a) those with odd numbers of external retarded indices, and b) the others with even numbers of external 
retarded indices. The n-point HTL vertex functions with one retarded index, which obviously belong to the first 
group a), are nothing but the HTL vertex functions that appear in the ITF, and vise versa Q. All the HTL vertex 
functions belonging to the first group a) are of 0{g^T^) , and satisfy among them the simple QED-type Ward- 
Takahashi identities, as in the ITF jl| |^. All the HTL vertex functions with a fermion pair belong to this first group 

a) . Those vertex functions belonging to the second group b) never appear in the ITF, namely their existence is 
characteristic of the RTF, and their HTL's have the high temperature behavior of 0{g^T^), one-power of T higher 
than usual. Despite of this difference we could verify that those HTL vertex functions belonging to the second group 

b) also satisfy among themselves the QED-type Ward-Takahashi identities, thus guaranteeing the gauge invariance of 
the HTL's in the real time thermal QCD/QED. The group b) HTL's consist of the n-gluon vertex functions. 

Comments and discussion are in order. 

i) With the present results we can feel at ease to perform in the framework of real time thermal field theories the 
nonperturbative analyses in the HTL approximation on dynamical phenomena such as the chiral phase transition. As 
an application we derived the HTL resummed DS equation for the retarded fermion self-energy function, enabling us 
to study in QCD/QED the chiral phase transition at finite temperature, which is now under investigation pO| . 

ii) It is well kn own th at all the HTL co ntributi on in the ITF are of order 0{g^T^). Are there anything wrong in 



our results, Eqs. ( 111.8c ), ( III. 18c -g), and ( [II. 24c -g), contradicting to the previous results? The answer is NO. It is 
shown Q that the n-point functions in the ITF is nothing but the special combination of those in the RTF, namely, 
just the n-point functions in the real time physical representation with only one external retarded index, all the other 
(n-1) external indices being the advanced ones. This means that the n-point functions with more than two external 
retarded indices never have their counterparts in the ITF. Furthermore we have already had some results indicating 
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the existence of the order 0{g^T^) contributions in the RTF |||, |T^. Our results as above have clearly shown that 
such order 0{g^T^) contributions really exist in the real time physical representation. Of course these unusual HTL's 
appear only in the n-gluon vertex functions in QCD. In QED there are no 3-photon vertex functions and no 4-photon 
HTL vertex functions, as already shown. 



iii) Another point to be noted is that the results, say, Eqs. ( 111.8a -c) have nothing in trouble with the number of 
inde pende nt components which is in the present case only one. As we can easily see the HTL contributions satisfy 
Eq.([i!L2^, 



(m^';(p,Q))*. 



Remembering that Eq.(III.21:) holds among the full two point vertex functions and here we are studying the HTL 
contributions, we can see that the additional power of T in Eq.([H.6c) comes from the boson distribution function, 
guaranteeing the number of independent components in the present case being only one. The same is true for 3- and 
4-point vertex functions with more th an two external retarded indices. 



iv) A s not iced a t the end of section IV, although the Ward-Takahashi identities ( IV. 21 ) and ( IV.Bg ), and also those 
( IV.Scf) and ( IV. 51 ) seem t o have a bit different structure from others in Eqs. (IV. 2), (IV. 3) and (IV. 5), all the identities 
( IV.2a| -g), ( IV.3a -g) and ( IV.5a -g) actually have the same structure. This fact can be clearly seen by noticing that, 
e.g., the right-hand-sides of Eqs. (IV.Sc) and (IV.Sl) are the HTL contributions to the {i?i?}-component of gluon 
polarization tensor which is actually the difference of two polarization tensors. 



(1 



2ns(go))sgn(qo)(n^'^ 



P + Q^O, 



namely, in the HTL approximation (where the external momentum —P = Q must be soft) 



su>^^'j,iP,Qy- 



2T 



RA 



{P,Q) 



<5H7^(F,Q)). 



(VI.l) 



(VI.2) 



Eq.( VI.2 ) shows that the right-hand-sides of Eqs. ( [V.5e| ) and ( [V.51 ) are nothing but the difference of two HTL gluon 
polarizatio n ten sors, t hus ha ving exactly the same structure as others. The same is true for the right-hand-sides of 
identities (IV. 21) and ( IV.3gD , though less clear but can be easily seen with explicit manipulations. 

v) We showed explicitly that in the physical representation in the RTF there exist those vertex functions with two 
external retarded indices having the high temperature behavior of 0{g^T^). This should be true for any n-gluon 
vertex functions with even number external retarded indices. 
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^ -i5£(P,Q) 



P+Q=0 

FIG. 1: One- loop diagram for the fermion self-energy and the definition of its HTL, JE. = denotes the equality in the HTL 
approximation. 




FIG. 2: One- loop diagram for the gauge boson polarization tensor andthe definition of its HTL, 511''''. In QED only the 
diagram (a) exists. 




FIG. 4: One-loop diagrams for the fermion-gauge boson vertex function. In QED only the diagram (a) exists. 




FIG. 3: Definition of the fermion-gauge boson HTL resummed vertex function *r''. 
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P,a,|Li Q,b,v 



FIG. 5: Definition of the 3-gauge boson HTL resummed vertex function . 




FIG. 6: One-loop diagrams for the 3-gauge boson vertex function. Figure (d) represents all possible diagrams with independent 
configuration with respect to the three external legs with momenta P, Q and R. 



S,\i R,v 




{g^Cpl*r^V(pQRS); Jqcd 1 

P+Q+R+S=0 



FIG. 7: Definition of the fermion-pair-2-gauge boson HTL resummed vertex function *r''''. 



S R S R S Q 




FIG. 8: One-loop diagrams for the fermion-pair-2-gauge boson vertex function. Their possible exchanged diagrams between 
the external gauge boson legs with momenta R and S should also be added. In QED only the diagram (a) and its exchanged 
diagram exist. 
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S,c,a R,c,p 




-ig'5''N:r^''''''(P,Q,R,S); 

P+Q+R+S=0 



P,a,^ Q,b,|a 

FIG. 9: Definition of the 4-gluon HTL resummed vertex function *Y'^^'^f"^ . Color indices of the external gluons with momenta 
R and S is summed over, which is indicated by the same color index c. 



S R S R S R 




FIG. 10; One-loop diagrams for the 4-gluon vertex function. For figures (a), (b) and (c), their possible exchanged diagrams 
among the external legs with momenta Q, R and S should also be added. Figures (d) and (c) represent all possible diagrams 
with independent configuration with respect to the four external legs with momenta P, Q, R and S. 




P Q 

FIG. 11: Possible contribution from the fermion-pair-2-gauge boson HTL vertex function to the HTL resummed DS equation 
for the fermion self-energy function Hr. denotes the corresponding vertex function and propagator being the HTL resummed 
ones. 



